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Abstract — In  this  paper,  we  discuss  the  development  and 
experimental  implementation  of  a  nonlinear  control  design  for 
magnetostrictive  transducers  operating  in  hysteretic  regimes. 
The  hysteresis  and  constitutive  nonlinearities  are  characterized 
using  a  homogenized  energy  framework  based  on  energy 
relations  at  the  lattice  level  employed  in  combination  with 
stochastic  homogenization  techniques  that  incorporate  material 
and  held  nonhomogeneities.  Using  this  framework,  we  employ 
nonlinear  optimal  control  theory  to  construct  open  loop  inputs 
for  tracking.  We  subsequently  employ  Pi-based  perturbation 
feedback  to  ensure  robustness  with  respect  to  model  uncer¬ 
tainty  and  sensor  noise.  Experimental  implementation  results 
at  frequencies  up  to  1000  Hz  demonstrate  the  feasibility  of 
the  method  for  high  speed  tracking  while  operating  in  highly 
nonlinear  operating  regimes. 

I.  Introduction 

Present  and  emerging  automotive,  industrial,  aeronautic, 
aerospace  and  biomedical  applications  require  actuators  that 
provide  large  force,  high  accuracy,  moderate  stroke,  broad¬ 
band  capabilities.  Moreover,  the  actuators  often  must  be 
compact,  reconfigurable  and  multifunctional.  For  a  num¬ 
ber  of  applications,  transducers  employing  ferroelectric  or 
ferromagnetic  drive  elements  meet  these  criteria  and  are 
being  considered  for  high  performance  control  systems.  Two 
common  ferroelectric  and  ferromagnetic  actuator  materials 
are  lead  zirconate  titanate  (PZT)  and  Terfenol-D. 

Whereas  these  compounds  provide  unique  actuator  capa¬ 
bilities,  they  also  exhibit  hysteresis  and  constitutive  nonlin¬ 
earities  that  must  be  accommodated  in  models,  transducer 
designs,  and  model-based  control  designs  to  achieve  strin¬ 
gent  tracking  and  control  specifications.  The  fact  that  the 
hysteresis  and  nonlinearities  are  rate,  stress,  and  temperature- 
dependent  compound  the  challenge  of  high  performance 
actuator  design. 

The  most  obvious  way  to  limit  hysteresis  and  nonlinear 
effects  is  to  restrict  devices  to  low  or  moderate  drive  regimes. 
Whereas  this  is  feasible  for  some  applications,  it  excludes 
the  devices  from  many  high  performance  applications  where 
can  prove  advantageous  over  traditional  actuators.  Secondly, 
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feedback  can  be  used  to  linearize  the  material  behavior  for 
certain  operating  regimes.  However,  it  is  illustrated  in  [9]  and 
Section  IV  that  the  authority  of  feedback  algorithms  degrades 
at  high  frequencies.  PZT-based  devices  can  also  be  linearized 
to  a  certain  degree  through  the  use  of  current  or  charge- 
controlled  amplifiers  [6],  [7],  However,  this  can  prove  expen¬ 
sive  when  compared  with  more  traditional  voltage-controlled 
amplifiers,  and  current-control  is  ineffective  for  maintaining 
a  fixed  dc  bias  as  required  for  various  applications  —  e.g., 
maintaining  a  fixed  position  with  an  x-stage  while  sweeping 
with  a  y-stage  in  an  atomic  force  microscope. 

In  this  paper,  we  experimentally  demonstrate  the  feasi¬ 
bility  of  a  nonlinear  control  design  for  high  speed,  high 
accuracy  tracking.  For  brevity,  we  focus  on  a  magnetic 
device  but  the  models  and  model-based  control  design  are 
equally  applicable  to  ferroelectric  materials  —  e.g.,  PZT- 
based  devices.  The  tracking  application  is  motivated  by 
recent  investigations  focused  on  the  use  of  Terfenol-D  trans¬ 
ducers  to  mill  automotive  components  at  high  speed  while 
maintaining  micron-level  tolerances;  however,  the  problem 
and  general  formulation  are  ubiquitous  for  a  wide  range  of 
applications. 

There  are  essentially  two  strategies  to  control  hysteretic 
actuators  using  nonlinear  models.  The  first  is  to  use  the  char¬ 
acterization  framework  to  construct  an  approximate  nonlinear 
inverse  that  linearizes  the  actuator  response  in  the  manner 
depicted  in  Figure  1(a)  —  e.g.,  see  [24],  Linear  control 
algorithms  are  then  used  to  achieve  control  objectives.  This 
method  has  the  advantage  of  linear  control  designs  but  has 
the  disadvantage  that  the  model  inverse  may  be  difficult  to 


(a) 


(b) 

Fig.  1.  (a)  Linear  control  design  employing  a  nonlinear,  model-based 

inversed  compensator,  and  (b)  nonlinear  model-based  control  design. 
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initiate  and  implement  at  high  speed.  Moreover,  scaling  or 
mapping  may  be  required  to  convert  amplifier  outputs  (e.g., 
voltage,  current,  or  field)  to  inverse  inputs  (e.g.,  magnetiza¬ 
tion).  Within  the  context  of  ferroelectric  and  ferromagnetic 
actuators,  this  approach  has  been  employed  using  Preisach 
models,  domain  wall  models  and  homogenized  energy  mod¬ 
els.  We  note  that  experimental  implementation  of  open- 
loop  and  feedback  designs  employing  Preisach-based  inverse 
compensators  are  reported  in  [3],  [8],  [2 1  ]— [23]  whereas  open 
loop  experimental  control  implementation  of  a  homogenized 
energy  inverse  are  reported  in  [4], 

The  second  strategy  is  to  construct  nonlinear  control 
designs  which  yield  input  signals  that  directly  incorporate 
actuator  nonlinearities  as  depicted  in  Figure  1(b)  —  e.g., 
see  [13],  [25].  If  optimal  control  theory  is  employed,  this 
approach  requires  the  approximate  solution  of  a  two-point 
boundary  value  problem  which  can  be  computationally  in¬ 
tensive.  However,  this  technique  avoids  the  real-time  imple¬ 
mentation  of  the  model  inverse  and  scaling  issues  associated 
with  inverse  compensators. 

In  this  paper,  we  demonstrate  the  experimental  imple¬ 
mentation  of  a  nonlinear  control  based  on  the  homogenized 
energy  framework  for  characterizing  hysteresis  in  ferro¬ 
magnetic  materials.  Open  loop  control  inputs  are  computed 
offline  through  the  approximation  of  a  two-point  boundary 
value  problem  derived  using  optimal  control  theory.  To 
provide  robustness  with  regard  to  model  and  measurement 
errors,  one  can  either  linearize  about  the  optimal  state  and 
control  to  construct  a  linear  perturbation  model  that  is 
amenable  to  LQR  design,  or  employ  direct  PI  feedback  on 
measured  perturbations  about  the  optimal  state.  We  employ 
the  latter  technique  since  implementation  in  this  manner 
provides  the  efficiency  of  classical  PI  designs  while  ac¬ 
commodating  frequency-dependent  hysteresis  and  nonlinear 
material  dynamics  via  the  nonlinear  open-loop  signal. 

The  hysteresis  characterization  framework  and  transducer 
model  are  summarized  in  Section  II.  The  nonlinear  open  loop 
and  perturbation  designs  are  discussed  in  Section  III,  and  the 
experimental  performance  of  the  technique  at  frequencies  up 
to  1000  Hz  is  illustrated  in  Section  IV. 

II.  Nonlinear  Constitutive  Relations  and 
Transducer  Model 

We  summarize  here  the  ferromagnetic  hysteresis  frame¬ 
work  developed  in  [16],  [20]  and  use  it  to  construct  a 
model  that  characterizes  the  dynamics  of  the  experimental 
magnetic  transducer.  As  detailed  in  [14],  [19],  analogous 
relations  quantify  the  hysteresis  inherent  to  PZT  and  hence 
the  framework  is  applicable  to  a  range  of  ferroic  compounds. 
It  is  noted  in  [4],  [14],  [18]  that  whereas  the  homogenized 
energy  framework  provides  an  energy  basis  for  certain  ex¬ 
tended  Preisach  models,  it  fundamentally  differs  from  clas¬ 
sical  Preisach  models  in  a  number  of  aspects  including  the 
direct  incorporation  of  thermal  relaxation  mechanisms  and 
the  capability  that  it  provides  to  characterize  noncongruent 
behavior. 


A.  Constitutive  Relations 

Letting  H  and  M  denote  the  magnetic  field  and  magneti¬ 
zation,  we  consider  the  lattice-level  Gibbs  energy  relation 


G(H,  M)  =  ip(M)  -  HM 


(1) 


where  the  Helmholtz  energy  is  given  by 

f  \rj{M  +  MRf  .  M  <  -Mr 

ip{M)  -  i  ~  mR )2  ) M  >  Mi 

[  HM!  -  Mr)  (C  -  Mr) 

(2) 

Here  M [ ,  M p>  and  ?/  respectively  denote  the  positive  in¬ 
flection  point,  the  local  remanent  magnetization,  and  the 
reciprocal  slope  after  switching. 

As  detailed  in  [14],  [16],  [20],  the  balance  of  G  and  the 
relative  thermal  energy  k.T /V,  where  k,  T  and  V  respectively 
denote  Boltzmann’s  constant,  temperature  in  degrees  Kelvin, 
and  a  reference  volume,  yields  the  kernel  relation 


M(H)  =  Z+  (M+)  +  x _  (M_)  (3) 


which  characterizes  hysteresis  at  the  lattice  level.  The  frac¬ 
tions  ,i;.|  and  x  of  positively  and  negatively  oriented  dipoles 
are  quantified  by  the  differential  equations 


x+  =  —p+-X. |_  +P-+ X- 
X-  =  —p-  +  X-  +P+-X+ 


(4) 


and  the  expected  magnetizations  due  to  positively  and  neg¬ 
atively  oriented  moments  are 


<M+) 

(M_) 


Q  Me~G^H’M^kTdM 
Q  e-Gtw^dM  ’ 

f~Ml  Me~G(H’M^kTdM 

Mi  e_G(HiM)v/kT dM 

J  — OO 


(5) 


The  likelihoods  of  switching  from  positive  to  negative,  and 
conversely,  are  given  by 


1  d-,e~G{H'M)V'kTdM 

P+ ~  =  TJt)  f^i_e  e~G^H'M^v !kT dM  ’ 

1  e~G^H'M)v'kTdM 

p~+  ~  T(T)  r-Mi+e  -G(H,M)V/kTdM 

J  —  OO 


(6) 


where  e  is  taken  to  be  a  small  positive  constant.  The 
relaxation  time  T  is  the  reciprocal  of  the  frequency  at  which 
dipoles  attempt  to  switch  —  see  [14],  [16],  [20]  for  details. 

To  incorporate  the  effects  of  polycrystallinity,  material 
nonhomogeneities,  and  variable  interaction  fields,  we  assume 
that  local  coercive  and  interaction  fields  are  manifestations  of 
distributions  having  associated  densities  vi(Hc)  and  v2(Hj). 
This  yields  the  macroscopic  field-magnetization  relation 


M(H) 


M(H+Hr,  tfc)ui(iTc)u2(fT/)dffcdff/. 

(7) 


For  the  device  characterization  results  summarized  in  Sec¬ 
tion  IV,  values  for  the  discretized  densities  were  estimated 
using  the  techniques  detailed  in  [14],  [17], 

B.  Transducer  Model 

To  model  the  experimental  transducer  depicted  in  Figure  2, 
we  assume  that  the  Terfenol-D  rod  has  area  A,  length  £, 
Young’s  modulus  Y,  density  p,  and  Kelvin- Voigt  damping 
parameter  C.  The  longitudinal  displacement  of  the  rod  tip  is 
denoted  by  Uf  . 

It  is  shown  in  [14]  that  the  uniaxial  stress  a  is  related  to 
the  magnetization,  given  by  (7),  via  the  constitutive  relation 

a  =  Y£  +  Ci-a1(M-Mo)-a2(M-M0)2  (8) 


The  model  can  subsequently  be  written  as  the  first-order 
system 

x(t)  =  Ax(t)  +  B(H(t)) 

(12) 

z(0)  =  x0 


where  x(t)  =  [ue(t),  iie(t))T,  xq  =  [uq,ui)t  and 

A  = 


0  1 
—k/m  —c/m 


B(H)  =  \ai(M(H)  -  M0)  +  a2(M(H)  -  M0f 


1/m 


In  the  subsequent  control  formulation,  the  input  is  taken  to 
be  u  =  H  or  u  =  nl  where  /  denotes  the  current  applied  to 
the  solenoid  and  n  is  the  number  of  coils  per  unit  length. 


where  Y  is  the  elastic  modulus  at  constant  magnetization,  £ 
is  the  longitudinal  strain,  C  is  the  Kelvin-Voigt  damping 
parameter,  a\  is  the  piezomagnetic  coefficient,  and  a2  is 
the  magnetostrictive  coefficient.  The  bias  magnetization  Mq 
includes  the  effects  of  a  permanent  magnet  and  the  initial 
magnetized  state  of  the  material. 

We  make  the  assumption  that  strains  are  uniform  and  given 
by 

*>  =  U~f-  (9) 

Balancing  the  forces  oA  for  the  rod  with  those  of  the 
restoring  mechanism  yields  the  lumped  model 


.  „d2ut .  .  CAdue,  .  YA  , 

d2Uf  .  .  due,  ,  .  ,  .  , 

=  -me-; -n-{t)  -  ce—{t)  -  kue(t) 

dtz  dt 

+Aai[M(H(t))  -  Mo)  +  Aa2[M(H(t))  -  M0] 


or,  equivalently, 

+  c~ 

=  ai[M(H(t))  -  M0]  +d2[M(H(t))  -  M0 ]2 


where 


III.  Perturbation  Control  Formulation 


A.  Open  Loop  Control 

We  summarize  here  the  open  loop  optimal  control  for¬ 
mulation  detailed  in  [10],  [12],  We  let  r  denote  a  reference 
signal  to  be  tracked,  y(t)  =  Cx(t)  denote  observations  and 
let  e(t)  =  Cx(t)  —  r(t )  designate  the  error.  The  augmented 
penalty  functional  is  taken  to  be 


J  =  \[Cx{tf)  -  r{tf)]T P[Cx(tf)  -  r{tf )] 

+  f  [H  —  XT  (t)x(t)]  dt 
Jt0 


(13) 


where  A  denotes  the  adjoint  variable,  the  Hamiltonian  is 


TL  =  i[e(£)TQe(f)  +  uT  (t)Ru(t)} 
+At  [. Ax(t )  +  [_B(it)](i)] , 


(14) 


and  Q,  R  respectively  penalize  large  errors  and  control 
inputs. 

Enforcement  of  necessary  conditions  to  minimize  (13)  — 
see  [2],  [5]  —  yields  the  control  input  relation 


along  with  the  two-point  boundary  value  problem 


to  =  pA£  +  me 


a\  =  Aa\  ,  a2  =  Aa2. 


(11) 


The  initial  conditions  are  ue( 0)  =  Uq  and  ^-(0)  =  U\.  The 
magnetization  M  is  specified  by  the  model  (7). 


Fig.  2.  Magnetic  transducer  design  employed  in  the  control  experiments. 


z{t)  =  F(t,z) 

where  z  =  [x,  A]T  and 

[  Ax(t)  +  [B(u)}(t) 

Fit ,  z)  = 

_  -ATX(t)  -  CTQCx{t )  +  CTQr(t) 


(16) 


(17) 


To  approximate  the  solution  to  (16),  we  employ  a  finite 
difference  discretization  defined  on  the  grid  tj  =  jAt,  where 
At  =  j?  and  j  =  0,  •  •  •  ,  TV.  Letting  Zj  ~  z(tj),  this  yields 
the  discrete  system 

[zj+ i  ~  Zj)  =  —  [f  (tj,  Zj)  +  F(tj+ 1,  Zj+ 1)] 

E0z0  =  [xo,  0]T 
EfzN  =  [0  ,-CTPr{tf))T. 


(18) 


The  solution  of  (18)  can  be  expressed  as  the  problem  of 
finding  Zh  =  [zq,  ■  ■  ■  ,  zn]  which  solves 

Hzh)  =  0.  (19) 

A  quasi-Newton  iteration  of  the  form 

4+1=*h+&,  (20) 

where  ^  solves 

E'(4)th  =  (21) 


is  then  used  to  approximate  the  solution  to  (19).  The  Jacobian 
has  the  form 


F{*h) 


Sq  R0 
Si  Ri 


where 


Eq 


Si  = 


1 

A t 


I  0 
0  I 


Sn-i  Rn-i 
Ef 

1 

2 


mB  k* 


~CTQC  - AT 


(22) 


■  (23) 


The  representation  for  Ri  is  similar. 

Remark  3.1:  It  is  shown  in  [13]  that  an  analytic  LU 
decomposition  can  be  determined  for  E(z^).  This  signifi¬ 
cantly  reduces  memory  requirements  and  is  fundamental  for 
efficient  solution  since  it  reduces  the  dimension  of  the  linear 
system  solution  to  8.  For  the  results  reported  in  Section  IV, 
solution  of  (19),  and  hence  computation  of  a  nonlinear 
control  input,  took  on  the  order  of  7  seconds. 


B.  Perturbation  Feedback 


The  relation  (15)  provides  an  open  loop  control  signal 
that  is  optimal  for  given  Q,R  in  the  absence  of  model  or 
measurement  error.  To  provide  robustness  with  regard  to  such 
uncertainties,  we  consider  perturbation  feedback  as  detailed 
in  [2],  [5],  Linearization  about  the  optimality  system  yields 

8x(t)  =  ASx(t)  +  B5u{t) 

Sy(t)  =  C8x(t)  (24) 


Sx(0)  =  xo 

where  5u,6x  and  5y  are  first-order  variations  about  u*,x* 
and  y*. 

To  facilitate  experimental  implementation,  we  use  classical 
PI  control  to  compute  Su;  that  is,  we  take 

du(t)  =  —Kpe(t)  —  Ki  f  e(s)ds.  (25) 

Jo 

The  final  control  input  is  then 

u(t)  =  u*(t )  +  Su(t).  (26) 

Remark  3.2:  The  computation  of  the  open  loop  input  u*  is 
performed  off-line  and  hence  does  not  affect  implementation 
speed.  The  computation  of  8u  can  be  performed  at  the  same 
speed  as  classical  PI  implementation.  Hence  the  experimental 
implementation  of  the  perturbation  control  is  as  efficient  as 
standard  PI  implementation. 


IV.  Experimental  Results 

The  control  algorithm  was  implemented  on  a  magne- 
tostrictive  transducer  having  a  Terfenol-D  rod  of  length  0. 1 
m  and  diameter  of  0.0125  m.  The  preload  on  the  rod  was 
10-14  MPa  and  it  was  subjected  to  a  magnetic  field  bias  of 
approximately  40  kA/m.  Strains  at  the  end  of  the  rod  were 
measured  with  a  capacitive  sensor  having  a  sensitivity  of  2.5 
mV  and  a  bandwidth  of  12.5  kHz. 

To  ensure  that  the  transducer  was  initialized  to  the  max¬ 
imal  remanent  value  prior  to  characterization  and  control 
experiments,  each  data  set  was  initiated  by  a  half-cycle  of 
a  1  Hz  sine  wave  having  an  amplitude  of  1  volt  (this  cor¬ 
responds  to  a  current  of  4.6  amp).  Sinusoidal  input  voltages 
and  corresponding  output  displacements  were  subsequently 
collected  at  100,  200,  300  and  500  Hz.  As  illustrated  in 
Figure  3,  the  degree  of  hysteresis  exhibited  by  the  device  is 
highly  frequency-dependent  thus  necessitating  the  incorpora¬ 
tion  of  magnetic  moment  dynamics  via  the  relations  (3)  and 
(4),  and  the  rod  dynamics  modeled  by  (10).  Details  regarding 
the  characterization  of  the  nonlinear  voltage-current  relation 
are  provided  in  [9], 

To  construct  the  model,  the  parameters  Mr,r),T,  7, 
ai,a2,m,c,k  and  densities  and  02  were  estimated 
through  a  least  squares  fit  to  the  combined  set  of  100,  200, 
300  and  500  Hz  data.  The  resulting  model  fit  is  plotted 
in  Figure  4.  The  model  with  this  parameter  set  was  then 
employed  in  all  control  experiments. 

To  illustrate  the  tracking  capability  of  the  nonlinear  control 
design,  experiments  were  performed  at  100,  200,  300,  500, 
700  and  1000  Hz.  The  amplitude  of  a  sinusoidal  reference 
signal  was  chosen  to  be  300  ppm  which  represents  an  operat¬ 
ing  regime  in  which  hysteresis  and  constitutive  nonlinearities 
are  significant  —  see  Figure  3. 

The  open  loop  and  perturbation  feedback  control  results  at 
200  Hz  are  respectively  plotted  in  Figures  5  and  6  whereas 
analogous  results  at  1000  Hz  are  plotted  in  Figures  7  and  8. 
For  both  frequencies,  it  is  noted  that  whereas  the  open  loop 
control  (15)  provides  reasonable  accuracy,  the  inclusion  of 
the  perturbation  feedback  term  (26)  improves  the  accuracy 


Fig.  3.  Frequency-dependent  current- strain  data  at  100,  200,  300  and 
500  Hz. 


(a) 


(b) 

Fig.  4.  Model  fit  to  the  experimental  data  at  (a)  200  and  (b)  500  Hz. 

of  both  the  magnitude  and  phase.  We  point  out  that  the 
characterization  experiments  did  not  include  1000  Hz  data; 
hence  the  accuracy  of  the  nonlinear  model-based  control 
at  this  frequency  illustrates  the  predictive  capability  of  the 
model. 

To  provide  an  initial  comparison  between  the  nonlinear 
control  technique  discussed  here  and  classical  PI  designs, 
the  tracking  authority  achieved  with  PI  inputs  at  1000  Hz 
is  illustrated  in  Figure  9.  At  this  frequency,  neither  the 
correct  amplitude  nor  phase  are  achieved  thus  yielding  errors 
on  the  order  of  the  reference  signal.  It  is  detailed  in  [9], 


Fig.  6.  Tracking  authority  provided  by  the  perturbation  feedback  control 
(26)  at  200  Hz. 


where  a  complete  comparison  between  the  two  designs  at 
multiple  frequencies  are  reported,  that  with  the  10  kHz 
sample  frequency,  PI  control  is  viable  at  frequencies  below 
approximately  500  Hz.  At  higher  frequencies,  the  nonlinear 
control  framework  proves  advantageous  since  it  provides 
high  accuracy  while  retaining  the  implementational  effi¬ 
ciency  of  classical  PI  designs. 


Fig.  7.  Nonlinear  open  loop  tracking  authority  provided  by  (15)  at  1000  Hz. 


Fig.  8.  Tracking  authority  provided  by  the  perturbation  feedback  control 
(26)  at  1000  Hz. 


Fig.  5.  Nonlinear  open  loop  tracking  authority  provided  by  (15)  at  200  Hz. 


Time  (s) 


Fig.  9.  Tracking  authority  provided  by  PI  control  at  1000  Hz. 


V.  Concluding  Remarks 

In  this  paper,  we  report  initial  experimental  results  demon¬ 
strating  the  tracking  authority  provided  by  a  nonlinear  control 
design  for  magnetostrictive  transducers  operating  in  highly 
nonlinear  and  hysteretic  regimes.  The  control  design  is 
comprised  of  a  nonlinear  open  loop  component,  based  on  a 
homogenized  energy  model,  employed  in  combination  with  a 
perturbation  PI  feedback  component.  Because  the  open  loop 
gains  are  computed  off-line,  the  closed  loop  implementation 
provides  the  efficiency  of  classical  PI  control  but  with  signifi¬ 
cantly  enhanced  tracking  capability  at  high  frequencies  since 
the  model-based  open  loop  control  incorporates  frequency- 
dependent  hysteresis  and  constitutive  nonlinearities.  This 
significantly  extends  the  results  of  [10],  [12]  which  focused 
solely  on  the  numerical  implementation  of  the  control  design. 

The  present  experimental  results  illustrate  sinusoidal  track¬ 
ing  at  frequencies  up  to  1  kHz  for  a  magnetic  actuator.  How¬ 
ever,  the  unified  nature  of  the  characterization  framework 
[14],  [19]  permits  direct  extension  of  the  theory  to  ferro¬ 
electric  (e.g.,  PZT)  and  ferroelastic  (e.g.,  SMA)  devices.  The 
numerical  implementation  of  an  analogous  control  design  for 
PZT  stages  is  illustrated  in  [11]  for  a  variety  of  reference 
signals  and  the  experimental  implementation  of  PZT  control 
design  is  under  present  investigation. 
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